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The spiral magnetic order in cubic MnSi-type crystals is considered using the model of classical 
Heisenberg ferromagnetics with an extra interaction of the Dzyaloshinskii-Moriya (DM) type be- 
tween neighboring atoms. For all shortest Mn-Mn bonds the DM vectors D are expressed via the 
D-vector for an arbitrary chosen bond, in accordance with the P2i3 space group of these crystals 
and bond directions. The D-vectors have the same length and 24 different directions so that the 
components of the D-vectors are changing similarly to the components of the corresponding bond 
vectors. After minimization of the magnetic energy it is found how the wave vector k of magnetic 
helices depends on D: for any helix direction k — 2(D X — 2D y — D Z )/(3J) where J is the exchange 
interaction constant and fc-vector is in the units of the lattice constant. The wave number k deter- 
mines both the sign and strength of global spiraling whereas locally, within a unit cell, the helical 
order can be strongly frustrated so that the twist angles between neighboring ferromagnetic layers 
may be even of different signs. Conical deformations of helices caused by an arbitrary directed ex- 
ternal magnetic field is also considered within the same model. The critical field of helix unwinding 
is found and it is shown that even in the unwound state there remains a residual periodic splay 
of magnetic moments with the splay angle proportional to (D x + D z )/J which can be measured 
by diffraction methods. The third component of the D-vector, D x + D y — D z , appears only in 
higher approximations. It is also demonstrated how the usually used continuous picture of moment 
distribution can be obtained from the discrete one in a coarse grain approximation. 

PACS numbers: 75.25.-j, 75.10.Hk 



I. INTRODUCTION 

The spiral arrangements appear in very different systems of condense matter physics and biology [l[ providing 
an important example of a primitive self-organization. Sometimes this self-organization results in rather intricate 
structures (may be the first staircases to more complicated biological patterns) like the blue phases in chiral liquid 
crystals or the Skyrmion textures in chiral magnetics with B20 crystal structure (MnSi, MnGe, Fei_j.Coj.Si) 

[Jf]]. There is a remarkable similarity in the physics of those so different systems [3]. Therefore it is interesting 
and challenging to look for local interactions responsible for macroscopic spiral ordering and, in particular, to find 
relationships between parameters of the local interaction and the final structures. For instance, the first question is 
what determines the period and handedness of helices? 

Most frequently, the helical arrangements are described phenomenologically, just by adding chiral terms like D n • 
curl n to the free energy of the system where vector n(r) describes a local system anisotropy. The period and 
handedness of helices are determined by the absolute value and sign of the pseudoscalar coefficient T>; of course this 
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term should be allowed by the system symmetry (lack of inversion is required). In the chiral liquid crystals, the value 
of T> is believed to be determined by steric interactions between neighboring chiral molecules but any quantitative 
theory is still absent because of baffling complexity of the problem. In cubic magnetic crystals of the MnSi type, the 
helical ordering is induced by the Dzyaloshinskii-Moriya (DM) interaction between neighboring magnetic moments Sj 

1CH13| can exist both in centrosymmetric and non-centrosymmetric crystals but only 



aa 



and Sj Is [9| . This interaction 



in the latter case it results in long-range helical structures. The most simple form of the DM term in the interaction 
energy is written as • [s, x sj] where pseudovector coupling coefficients D.^ originate from superexchange interaction 
via neighboring atoms. 

It was shown by Moriya that vectors are determined by asymmetry of local atomic arrangements. Recently 
it was carefully investigated, using synchrotron radiation and polarized neutron small-angle diffraction, what is the 
relationship between the handedness of the MnSi crystal structure and the sign of magnetic helix [if]] . It was 
unequivocally demonstrated that the structural chirality rigorously determines the magnetic chirality of these com- 
pounds. The interplay between structural chirality is also important for spin chirality in multilayer structures 
However, no analytical formula related T> and has been found yet even for the simple case of interaction between 
magnetic moments of neighboring atoms. This problem is indeed even more complicated because the itinerant nature 
of magnetism in MnSi crystals. Detailed numerical studies of the DM interaction in MnSi-type crystals are performed 
by Hopkinson and Kee [18| and their results will be discussed below in comparison with the results of the present 
paper. 

A generic feature, intrinsic to complicated spacial structures, is frustration of their local ordering: the configurations 
with the best local energy cannot be spread for all the system. In the chiral liquid crystals and MnSi heliomagnetics, 
the most favorable double-twist local arrangements are restricted by topological constrains existing in two-dimensional 
and three-dimensional spaces 

aa 

. Nevertheless, the stable three-dimensional blue phases are possible in a narrow 
range of temperatures whereas the two-dimensional system of spin vertices can be stabilized by external magnetic 
field. However, in complicated crystals like La2Cu04, there is an additional mechanism of frustrations [jjlH] related 
with different orientations of Dy for different inter-spin bonds. We have found that this type of frustration is also 
intrinsic to the MnSi structure: at the atomic scale, spin rotation in the magnetic helices can be strongly nonmonotonic 
because of different orientations of for different bonds. As a result, unequal righthand and lefthand twists usually 
coexist in each of the helices. For instance, the coefficient T> can become zero even for nonzero Dy. 

In this paper we first find a simple analytical expression related the helix period and handedness with the DM 
vectors of the shortest Mn-Mn bonds. In this consideration, we adopt a usually used approximation that the exchange 
interaction is stronger than DM interaction and twist angles between neighboring atoms are correspondingly small. 
It is demonstrated that for general orientation of DM vectors the MnSi helical structure is frustrated so that the twist 
angles between neighboring ferromagnetic layers have different signs. Then, using the same approximations, conical 
distortions of helices in external magnetic fields are analyzed and the critical field of helix unwinding is determined. 
Finally, the parameters of the known phenomenological description of MnSi helices are related to the DM vectors of 
the microscopic theory. 
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II. MAGNETIC ENERGY AND CRYSTAL STRUCTURE 

In a simple model of classical Heisenberg ferromagnetics with an extra interaction of the Dzyaloshinskii-Moriya 
(DM) type, the energy of the system is expressed as a sum of pair interactions between magnetic atoms and the 
interaction of individual atoms with an external magnetic field H: 

E = - y^(-JijSi ■ Sj + T>ij ■ [si x Sj]) - gun H Sj, (1) 

where the first summation is taken over all pairs {ij} of magnetic atoms, Jy is the exchange interaction constant, 
T)ij are the Dzyaloshinskii-Moriya vectors, and g/is is an effective magnetic moment. Hereafter we will use the 
approximation of classical spins with ferromagnetic exchange coupling between nearest neighbors, Jy = J > 0, and 
suppose that | | = 1 for all magnetic atoms. Using the nearest neighbor approximation for the MnSi-type crystals 
we can rewrite Eq. ([I]) as 

1 6 

E = - ^2^2(-Jsi ■ Sj + Bij ■ [si x Sj]) - .g/i B ^H • m, (2) 

i j = l i 

where the first sum (i) is taken over all Mn atoms of the crystal and the second sum (j) is over six nearest Mn 
neighbors of i-th atom, i.e. over six bonds connecting each Mn atom with its nearest neighbors. 

It should be emphasized that vectors Dy have different orientations for different bonds and relations between 
different Dy can be found from symmetry reasons. In the unit cell of MnSi-type crystals, the magnetic atoms are at 



the 4a positions with coordinates (x, x, x), (1 — x, \ + x, i — x), (5 — x, 1 — x, \ + x), {\ + x, \ — x, 1 — x) 2l|- There 
are 24 crystallographically equivalent bonds between nearest neighbors, and DM vectors for all bonds can be obtained 
from DM vector of any bond using symmetry transformations of the point group 23 of MnSi crystal augmented with 
the inversion; the latter appears owing to the permutation of atoms: Dji = -Dy. The procedure is, in principle, very 
simple: we choose an arbitrary vector (D x ,D y , D z ) as DM vector of the "first" bond, say (—2 ), connecting 

Mn atom at ( ) with Mn atom at ( X i 2 I X 5 2 X ). It should be an arbitrary vector because no symmetry 
operation is associated with the bond. Then, for this cubic space group, any shortest bond between Mn atoms can be 
obtained from the "first" one using cyclic permutations (i.e. threefold rotations) and/or sign inversions (i.e. twofold 
rotations and inversions) of its vector components. The components of the corresponding DM vector are expressed 
through D x , D y and D z by exactly the same cyclic permutations and/or sign inversions. The list of all the bonds and 
their D-vectors is given in Table HI All have the same length: |Dy| = D. 

It was shown [3,0] that the DM term should be accompanied by two extra terms: 

'D^fa-sj) is, •D )(D U -s,i\ 



2 V 4J 2 J 



(3) 



These terms are of the same order of magnitude as the DM term because for small DM interaction, |Dy| <§; J, small 
angles between and Sj are of order |Dy|/J. In Eq. @, with the same accuracy we can put = const for a 
small domain containing all 12 different interatomic bonds and averaging of (Dy) a (Dy) / g over these 12 bonds gives 
^D 2 6 a p. Thus Eq. ([3]) can be written as 



E ' = \H^-^ ( 4 ) 
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This term is isotropic and can be included into exchange term of Eq. ([T]). We see that in cubic crystals the extra 
terms given by Eq. ((3]) are not important if \D%j\ <C J- 

Having a list of DM vectors one can calculate and minimize numerically the magnetic energies (JlJ and © for rather 
big sets of atoms with distinct moment configurations, including those similar to Skyrmions. However in this paper 
we will use only analytic calculations for detailed consideration of helical structures. 

III. HELIX STRUCTURE IN THE ABSENCE OF MAGNETIC FIELD 

Let us start with the simple case of free helices which provides a minimum of the magnetic energies ([T]) and ([2]) for 
H = 0. As it is known from numerous experiments, well below the Curie point the magnetic moments of Mn atoms 
form a simple helix. However, in the general case, the helices may be different for four different types of equivalent 
atomic positions in the unit cell (i.e. for positions with four different orientations of the local threefold axes) , and this 
is confirmed by numerical simulations. Then, the spin structure can be described by the superposition of four helices 

s f = n* os cos k • r + n* in sin k • r, (5) 

so that the condition |s*| = 1 is obviously obeyed. Here r is an atomic position, n' os and n* in are perpendicular 
unit vectors, k is a wave vector. Upper index t = 1,2,3,4 enumerates four atomic positions in the "first" unit cell: 
(x, x, x), (1 — x, i + x, i — x), (h — x, 1 — x, i + x), (i + x, \ — x, 1 — x), correspondingly. We express all the atomic 
positions in the units of the cubic unit cell so that both r and k vectors are dimensionless. Notice that all four spirals 
given by Eq. ([5]) have the same wave vector whereas in the Skyrmion-like structures there is a set of different wave 
vectors. 

It can be easily seen that the spins of sublattice t rotate in the plane perpendicular to vector 

Kot = [n* os x n* si J, (6) 

whose direction generally speaking can differ from that of wave vector k. 
Now, using Eq. ([5]) we can rewrite expression ^ for energy (at H = 0) as 

Eheiix = \ Hi X^=i{(-^ n '°s • tSos + D, 3 • [n*oi x n'oi]) cos k ■ ri cos k ■ r, 
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(7) 



where upper index t(i) designates "type" of i-th atom or its belonging to one of four sublattices, t(i) — 1 . . .4. The 
sines and cosines products in Eq. ([7]) can be expressed through sines and cosines of arguments k ■ (2r^ + ) and k • , 
with Yij being the distance between z-th and j-th atoms, = rj — r^. It should be noted here that dependence on 
atomic positions is present in these trigonometric functions only. After summation over all unit cells of the crystal all 
the functions of argument k • (2r^ + r^) vanish because k is not a reciprocal lattice vector, whereas the functions of 
argument k • Yij do not depend on the unit cells to which i-th and j-th atoms belong. Thus we can rewrite expression 
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for energy as 

E helix = ^ Eti Ei= l {(-^nj M ' n - + ■ K os x n*^]) cosk ■ r. (J 
+(-J <os ■ n!to + D 2j • [n j. os x n*^ ] ) sin k • r i6 
-(-JKtn ■ n 'oi + D^i ' [n l si „ x n*^]) sink • r y 

+ (-^ n k» ' n lm + D ij ' K„ X COsk • r y}> 

where N ce u is the number of unit cells in the crystal, and first sum is now taken over four magnetic atoms within one 
unit cell, so t(i) = i. 

Eq. (|S]) seems to be simple enough for possible minimization. It depends on 15 variables only: three angles 
characterizing orientation of triad {n* os , n' irl , n* ot } for each of four Mn sublattices, and three components of the wave 
vector k. The exchange interaction constant J, DM vectors D y - and distances between neighboring atoms are 
supposed to be invariable; we do not consider phenomena (e.g. magnetostriction) changing these parameters. Below, 
instead of numerical minimization, we will look for approximate analytical solutions relying on a small parameter. 

Indeed, it is known from the experimental data that the spin-orbit interaction is usually much smaller than the 
exchange one, so that \Dij\ <C J. In the MnSi crystal, strong exchange interaction leads to the ferromagnetic ordering, 
whereas small DM interaction results in small canting of neighboring magnetic moments. Small values of D / J and 
|k| (the latter is supposed to be of the order of D / J) allows us to perform approximate minimization of the energy. 
In particular, we can take advantage of the fact, that in this approximation the triads {n* os , n* in , n* of } have close 
orientations for four Mn sublattices. 

Let us suppose that all triads {n* os , n* in , n* ot } are close to triad {ni, 112,113}, with 113 being a vector, which can 
differ from rik = k/fc (fc = ±|k| is a wave number; the sign defines the chirality of the helix), and {ni,ri2} being 
an orthonormal basis of the plane perpendicular to 113, [ni x 112] = 113. It is useful to introduce vector 113 7^ rik for 
considering all partial helices on an equal footing and for possible future consideration of anisotropy effects. 

The deviation of triad {n* os , n* in , n* ot } from {ni, 112, 113} can be considered as a sum of two small rotations: first 
one on the axis 113 by the angle Lo t , and another one by the angle jt -L 113. Then, the vectors n* os , n* in , and n* ot can 
be represented as 

n Ls = (-R*ni) costJt + (i?*n 2 ) sintJt, 

n sm = -(-R'ni)sinw t + (i?*n 2 ) cos uj u (9) 

n rot = -R' n 3, 

with i?* = i?(7t) being the rotation matrix by the angle 7*4. The substitution of Eq. ^ into Eq. (JSJ gives 

E heUx = ^ Eti E%i{-J{(R^i) ■ (tfWnx) + (R l n 2 ) ■ (i?<«n 2 )} cos(w 4 - w ty) - k • r y ) 
-J{(i?'n 2 ) • (i^V) ~ (R l nx) • (i^n^jsm^ - u t(i) - k • r tf ) 
+T> VJ ■ {[(i? J rn) x (JfOIni)] + [(i?*n 2 ) x (i?*«n 2 )]} cos(w< - u t{j) - k • r«) 
+T> VJ ■ {[(i?*n 2 ) x (JfO0 ni )] - [(i? l m) x (i? t (j) n2 )]}sinK - w t(j) - k • r< 3 -)}. 

The seeming dependence on arbitrary vectors ni and xi2 can be eliminated with use of following two equations 



(ni) a (ni) /3 + (n 2 ) a (n 2 ) / 3 = S a0 - (n 3 ) a (n s )p, 
(n 2 ) a (ni) /3 - (ni) a (n 2 ),3 = -e a / 3j (n 3 ) 1 . 



(11a) 
(lib) 
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Then, 



Ehelix = T,t=i Yfj=i{- JR ii* R %Z (<W ~ (n 3 ) Al (n 3 ) l ,) cos(w, - u t(j) - k • r tJ ) 



+JK t afi Ra£'e^ CT {n 3 ) c7 sin(cj l - w t(i) - k • r y ) 



(12) 

+(D ii3 -)a£ Q /3A^ A1 i?A^ ) (V - (n 3 ) AI (n 3 ) t/ ) cos(wi - u t{j) - k • r. y ) 

-(Di^aSapxR^^R*^/ ] e^ v<T {^)a sin(^ - W t y) - k • Yij)}. 



The phases u> are presented in expression (|12p only as differences, which is in agreement with the fact, that they 
are defined up to a constant term associated with the origin choice. The matrix R("f) can be expressed through the 
coordinates of the vector 7: 



(R(i)) afi =cos 7 ^-^e a ^7. + i ^7a7^ 



(13) 



~ S a /3 ~ £ a /3crla + \ 7/3 ~ S a ^ 2 ), 

with 7 = |7|. Then, we can rewrite the energy up to second order terms on |k| ~ (uii — ^t(j)) ~ 7 ~ {P/ J)' 

E heUx =N cell (-12J + e 1 +e 2 ), (14) 



with 



1 4 6 fj 1 

i=l 7 = 1 J 



(15) 



1 4 6 fj 1 

As we can see, the energy breaks up into two contributions depending on uj and 7 and related to each other only by 
direction of the axis n 3 . 

The routine minimization of ei as a function of independent variables (w% —u>i), (ui3 — u\ ) and (W4 — u>i ) (performed 
with Wolfram Research Mathematica) gives 

W2-W1 = -l(k x + k z )(l-8x)-^f±((n 3 ) x + (n 3 ) z ), 

W3-W1 = -i(fc x + ^)(l-8a;)- a^((n3) a + (n3)„), (17) 
W4-wi = -3(^ + ^)(l-8^)-^7^((n 3 )y + (n 3 ),), 

and 

/ (D x + D z ) 2 3 ;2 D a -2D V -D z \ 
ei = Jl j2 + -fc -/c (n k -n 3 )J. (18) 

It is natural that t\ does not depend on x, because real atomic coordinates are not involved in energy definition (fTJ. 
From the physical point of view it is more interesting that e\ is also independent of orientation of wave vector k 
relative to coordinate axes. This fact is in a good accordance with experimental evidence of the possibility of easy 
reorientation of helix axis in MnSi by magnetic field, surface anchoring or other effects. 

The minimization of the expression (1181) on k gives us the following value for the wave number depending on DM 
vector D and rotation axis orientation relative to k 

im n 2fD x -2D y -D z \ 

fc(D,n k n 3 ) = - f (n k -n 3 ), (19) 



3 V J 
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and the part t\ of unit cell magnetic energy can be finally expressed as 



ei = J 



(Ar + g Z 
J 2 



1 (gs 

3 



2£> v - £> z ) 5 



J 2 



(n k • n 3 ) 



-(k- fc(D,ri^K3)) : 



(20) 



thereby the elasticity of the helix is § J. 



It is evident from Eq. (I20p that ei has minimum, when n 3 coincides with nk, but it should be noted that, when 
n 3 deviates from by a small angle 7 ~ (D/J), the energy obtains only a small positive addition of the order of 
J{D/J) A . Thus, there is a soft mode associated with deviation of rotation axis from helix direction. This deviation 
gives a cycloidal component to the helix. 

Eq. (|16p for e 2 contains vectors 7* t only within differences, so, in analogy with lj, the vectors jt are defined up to a 
constant vector perpendicular to n 3 . However, unlike the case of cj, where a common term added to all phases results 
only in redefining of the origin, the constant vector added to all vectors j t turns rotation axes for all four sublattices 
by the same angle. In fact, this turning has been already taken into account in the dependence of e\ on (nk ■ n 3 ). 
We should also remind that Eq. (fT6|) has been obtained in the assumption that 7 ~ (D/J). In order to avoid the 
ambiguity of definition of vectors 74, we can use additional condition J3t=i 7* = 0- 

The routine minimization of 62 on independent variables (7*2 — 71), (73 — 71) and (74 —71) (performed with Wolfram 
Research Mathematica) gives 

'72-71 = ^^{(n 3 ) x (n 3 ) z -(n 3 ) 2 -(n 3 ) 2 , 

(n 3 ) x {n 3 ) y + (113)3,(113)3, (n 3 ) x (n 3 ) z - (n 3 ) 2 - (n 3 ) 2 }, 

73- 71 = ^f-{(n 3 ) x (n 3 ) y -(n 3 )l-(n 3 )l 
(n 3 ) x {n 3 ) y - (n 3 ) 2 - (n 3 ) 2 , (n 3 ) x (n 3 ) 

74- 7i = ^r^{(n 3 ) x (n 3 ) y + (n 3 ) x {n 3 ) z , 
(n 3 ) y (n 3 ) z - (n 3 ) 2 - (n 3 ) 2 , (n 3 ) y (n 3 ) 



(n 3 ) y (n 3 ) z }, 
(n 3 ) 2 - (n 3 ) 2 }, 



(21) 



and 



C-2 



J 



(22) 



The full energy of the helix can be now written as 



E helix = N cM J (-12 - 2 {Dx+ j2 Dz)2 - 2D J Dz) \ n k ■ n 3 ) 2 + \{k fc(D, fi^)) 2 ) , (23) 

where the contribution —12 J per unit cell is associated with ferromagnetic ordering of spins (12 is the number of 
Mn-Mn bonds per unit cell), the term —2 ( D:c+ J D ^ , combining e 2 with a part of e 1; describes the energy of divergence 
between sublattices (the energy of canting) , and remaining terms represent the twist energy of the helix. 
The minimum of the energy (1231) is achieved, when n 3 = n k . Then, 

*~ = jw (-n - 2<^£ - + » ( » _ mr ) , (24) 



where 



, 2 ( D x - 2D V - D z . 

fe(D) = -f— f -)• (25) 
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It is also useful to write the full divergence vectors for sublattices using the formula 

ip = ujn 3 + 7. (26) 

then, from Eqs. (|T7J| and fl2T} follows 

$2-$x = ^^^(e 2 -e 1 ) + ^(l-8a;)(k.e 2 -k.e 1 )n3, 

= ^^^(e 3 -e 1 ) + ^(l-8x)(k-e 3 -k-e 1 )n 3 , (27) 
= ^^±^(e 4 -e 1 ) + ^(l-8x)(k-e4-k-e 1 )n 3 , 

or, using the additional condition X)t=i V^t — 0, 

r V3 Ar + £, %/3 

^ = -4 J e 4 + — (l-8x)(k-e t )n 3 , (28) 

where ei, e2, e 3 and e 4 are the unit vectors directed along 3-fold symmetry axes corresponding to four sublattices: 

ei = ^(1,1,1), 

02 = ^(-1,1.-1). 

e 3 = ^(-1,-1,1), 
[e 4 = ^(1,-1,-1). 

Finally, we can express vectors n cos , n S i„ and n rot for each of the sublattices: 

n* os = m + [iJ t x m] = m + ^^±^[e t x m] + ^(1 - 8x)(k • e t )n 2 , 

n* m - n 2 + [<& x n 2 ] = n 2 + ^^±^[e t x n 2 ] - ^(1 - 8x)(k • e t )m, (30) 
n* ot = n 3 + $ x n 3 ] = n 3 + ^^±i[e ( x n 3 ]. 

Thus, in the common case each sublattice has its own rotation axis different from as it shown in Fig. HJa). The 
divergence of the axes is defined by Eqs. (|50"]) . The substitution of Eqs. ([50)1 into Eq. ([5]) gives the spin of an individual 
Mn atom 



(29) 



Si = (m + f^[e t(l) X m] + 3^(1 - 8x)(k • e t(i) )n 2 ) cosk • r, 
+ (n 2 + ^^±^[e i(i) x n 2 ] - - 8x)(k ■ e t(i) )m) sink • r, 



(31) 



The question can arise: why parameter x appears in Eq. (|31[) . in spite of its absence in the initial Eq. ([5])? In fact, x 
is brought in by Eq. ([5]) depending on real atomic positions. Indeed, parameter x is implicitly involved with atomic 
coordinates r,-, and expressing r as 

r = p + r', (32) 

with p being a period of the crystal or the origin of a unit cell, and r' £ {(x, x, x), (1 — x, h + x, | — x), (i — x, 1 — 
x, o + x), (| +3;, I - s, 1 - a:)} being an atomic position within the cell, and using the condition k-r' <C 1, we can 
rewrite Eq. (|3Tj) as 

Si = ( ni + ^^±^[e i(4) x m] + {^(1 -8a;)(k-e t(i) ) + k-r; (i) }n 2 ) cosk ■ p, 
+ (n 2 + ^^±^[e t(j) xn 2 ]-{|(l- 8a:) (k • e t(i) ) + k • r^} m) sink • Pi , 

where parameter x vanishes in curly brackets after substitution of corresponding e and r'. 
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Fig. [2] shows spin helix structures calculated with the use of Eq. (|3Tj) for directions (001), (111) and (Oil) of spiral 
axis in the case of 113 = rik. The wave number |k| is chosen to be equal to 27r/40, which is close to the experimental 
value for MnSi crystal. The canting D-vector component is taken to be rather large, D x + D z = 0.4J, in order to 
emphasize nonmonotonic frustrated behavior of helices. In this case 

\Dx + D*\ „ (8*-l)|k| r . 

U 8 1 ' 

and the main contribution in a canting is due to the terms with [e t (j) x ni^]. 

The extended sinusoids with phase shift of J represent spin projections on directions ri! and n 2 , the low ones do 
spin components along the wave vector, which appear due to the canting. The coloured circles correspond to magnetic 
atoms in four sublattices of the crystal. Being all different in the common case, the sinusoids connected with different 
sublattices can coincide for certain directions of k reflecting the symmetry. For instance, in the case of k || (001) 
(Fig. [U(a)) the extended sinusoids coincide, which correspond to sublattices 1, 3 and 2, 4; but the circles of different 
colors do not coincide because the atoms belonging to different sublattices lie in different planes perpendicular to 
k. On the contrary, in the case of k || (111) (Fig. [2jb)) coincide both the extended sinusoids and related circles 
corresponding to sublattices 2, 3, 4, because their sites lie in the same atomic planes perpendicular to k. Both two 
cases, (a) and (b), due to their symmetry ((001) is the axis 2i and (111) is the axis 3 of the crystal) are characterized 
by only two extended sinusoids, which defines the view of the black saws representing the angle </> between spins in 
the successive ferromagnetic planes, tan0 = (s ■ ri2)/(s • ni). The symmetry is also reflected in the low sinusoids 
corresponding to the spin components parallel to k. 

When k || (011) (Fig. [2c)) there are three extended sinusoids, which is also reflected in the view of the cooresponding 
black saw (Fig.[2jd)). In this case, coincide the extended sinusoids corresponding to sublattices 2, 3 and the low ones 
corresponding to sublattices 1, 4. It may seem from Fig. (5Jc) that sublattices 1 and 4 have the same projections on k, 
but it is evident from Fig. Hl^d) that it is not so. Indeed, the atoms of types 1 and 4 lie in close but different planes; 
their proximity is due to the approximate equality 4x ~ \ ■ 

The "saws" in the figures show a frustration which is the intrinsic feature of the helices due to the spin canting. 
Owing to the frustration the spins in successive ferromagnetic planes can rotate inversely to twist direction of the 
helix. 



IV. HELIX IN MAGNETIC FIELD 



Let us return to the case of nonzero magnetic field H. Numerous experimental data, as well as the results obtained 
with the use of phenomenological theory, show that in this case the magnetic moments of Mn atoms form a conical 
helix. The peculiarity of such structure is that the field of magnetic moment posesses a constant component directed 
along H. Then, the spin structure can be described by the following expression 

s* = n' ot sin £ + n' os cos £ cos k ■ r + n' in cos £ sin k • r , (35) 

which satisfies the condition |s'| = 1. Here £ is an angle characterizing the ratio between constant and rotating 
components of the magnetic moments. The substitution of Eq. (|35|) in the expression ([2]) and summation over all the 
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unit cells of the crystal give 



EhelixM = ^ Ei=l Ei=l(- J Kot ■ "rot + D ij ' [ n rot X n rot}) sin & Sin& W) 

+ Eli Ej =1 {(- •/„ • n<« + D y • [4 os x n^]) cosk • r y - 
• n*^ + D. (J • [n 4 cos x n*^]) sink • r ti 



-(- ^< in • n* ( J i + D, 3 • [n 4 sm x n'^]) sink • r„ 

+(- J<in ■ n l ( in + B ij ' Km x cos k ' r «j> cos & cos &(j) 



(36) 



(37) 



-N ceU gHB Ei=l H • n rot sin &' 

Using Eq. © we can rewrite the energy up to second order terms on |k| ~ (w, — w t(j)) ~ 7 ~ (D/ J)'- 

E helix , H = %^ Eti E'=i{-2Jcos(6 - & y) ) 

+2 (|(a>i - w t (j-) - k • Vijf - (Dy • n 3 )(cJi - cj t(j) - k • r. y )) cos& cob&c?) 

+ (7i - 7tb)) 2 - Dy • (7i - 7t(j))) (cos&cos£ t(j) + 2 sin& sin&y))} 

-N ce ugix B (& ■ n 3 ) Ei=i sin 6- 
The last term has been obtained supposing that g^sH ~ D 2 /J, which allows us to use approximation n* ot w n 3 . 
The first term in curly brackets has zero order on (D / J) and achieves the minumum, when 

^1=^2=^3=^4 = £. (38) 

Then, we can rewrite the energy as 

E he ii x ,H = N ce u(~12J + ei cos 2 $ + e 2 (l + sin 2 - 4s// B (H • n 3 ) sinf), (39) 

where ex and e 2 are the same as in the previous section. The minimum of the energy (|39p is achieved, when H || n 3 
and n 3 = n^. Then, 

T? AT jf 10 o (^ + ^) 2 1 (Ac - 2fl„ - A,) 2 g^gg . \ 

Eheltx.H = N ceU J 1-12 - 2 — j| cos z £ - 4 — — — sm£ I , (40) 



or, after minimization on £, 



t (D x +D z ) 2 l{D x -2D y -D z f l2{ gi i B H) 2 



E heltx ,H - N cell J ( -12 - 2 - - ______ ) , mi , 



and 

Qin £ — 

(D x - 2D y - D z ) 



. 6Jgn B H 
sm £ = 77^ 7TT7- ( 42 ) 



The helix vanishes, when £ = §j which corresponds to magnetic field 

(A* - 
,. = — — 

QJg/j, 



Hc= (D x -2D y -D z r (43) 



(this field, unwinding the helix, is usually called H C 2 22j). 

Thus, in the approximation used here, the wave vector k is directed along the applied magnetic field. The spin 
structure for each sublattice obtains a permanent component along its rotation axis n rot , which is different from 
direction nt || H, see Fig. [ljb). The helices become conical, and as the magnetic field increases up to H c , the 
permanent components grow proportionally to the field value, Eq. (I42p . whereas the rotated components decrease 
correspondingly. It should be noted that if the field is less than H c , it does not affect the axes n rot and they have the 
same directions, defined by Eq. (f2"5)l . as in the absence of magnetic field. 
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COARSE GRAIN APPROXIMATION 



In the phenomenological theory, the energy of chiral magnetics is usually written as 



E= I (^+PM.^M]|-/, 



dr m dr„ 

where M(r) is the local magnetic moment changing smoothly in space. Below it will be shown how this equation and 
phenomenological coefficients J and T> arise from microscopic consideration based on Eq. ([2]). 

Let us suppose that the typical scale at which a considerable change of magnetic moment occures is much bigger 
than the size of the unit cell, and any two neighboring spins are rotated relative to each another by a small angle of 
order (D/J). Let Si and Sj be the neighboring spins and Sj — R(Cij)si, or 



sin Qj .-» 1 — cos Ci 7 
s 2 cos Cij + t [Qj x Si\ H -2 -CijiCij ■ s i), (45) 



where Cij = \Cij\- Then, 



Bi ■ Si « 1 - $/2 + (at ■ C^) 2 /2 = 1 - $ j± /2, (46) 

[Si X Sj] PS £y - Sj(Sj • (y) = Cij-L, (47) 

where is the part of vector Cij perpendicular to unit vector Sj. 

If spins vary slowly along the crystal, we can take a domain, which is much bigger than a unit cell, but small enough 
as compared with the scale at wich a considerable change of spin occures. The spins within this domain are almost 
unidirectional. Consequently, all the vectors Cij± nearly lie in a plane. Let us neglect this "nearly" and suppose that 
all Cij-L belong to the plane perpendicular to some mean spin s. In addition, staying in the domain, we can discard 
the parts of Cij parallel to s and suppose that Qj± = Cij- 

Another simplification can be achieved with taking into account that a spin does not change upon round a close 
chain of atoms. This means that the sum of small angles Cij corresponding to the bonds of the chain should be equal 
to zero. For instance, for the triangle formed by spins Si, S2 and S3 we obtain the following condition 

C12 + C23 + Csi - 0, (48) 

which is evidently satisfied, if Cij = Cj ~ Ci: where Ci and Cj are perpendicular to s vectors associated with i-th and 
j-th atoms. 

With the approximations accepted, using the expressions (|46|) and (j47| . the energy ([2]) in the absence of magnetic 
field can be now rewritten as 

E = \ E E(- J + 5(6 - ^ ) 2 + D « • (5 - 5))- (49) 

i j=l 

The minimization of the energy on the vector Q associated with i-th atom and belonging to the plane perpendicular 
to s gives 

dE 6 

^ = - E( J (0 - 5) + D ij± ) - 0, (50) 
dQ 
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1 6 1 ( 6 \ 1 6 

i=i \j=i / ± j=i 



- 1 P a + £>s 

\/5 ^ Bm± ' 



where symbol "_L" designates the projection of a vector to the plane perpendicular to s, etu) is one of the vectors 
(|2T)1) . The last therm in the right part of (I5T1) reflects the evident fact that the sum of DM vectors associated with all 
the bonds of an atom should be parallel to the 3-fold axis to which the atomic position belongs. 
Let us introduce a new vector Q connected with the old one by equation 

Q = t+Ae m ±- (52) 

where A is a constant. Then, 

Ci = 6 g c > + v 3 + 7! 1 J et «^' (53) 

where it has been taken into account that an atom belonging to one of four sublattices has two neighbors at each of 
three remaining ones, so 

6 

53e tW = -2e t(i) . (54) 



Choosing 



we can now write 



A =4^1, (55) 



7 ~ 7 , y/3 D x + D z 

Ci = Ci + —> -j e t(l)x , (56) 



S = j?£§> (57) 



3=1 

and, after the substitution of Eq. (|56|) into Eq. (|49|) . 



IE* E 5 =1 {- J + 1(0 - CO 2 + D« • (0 - 5)} 



! — ./ + — I I : — I I" -!- I I : : • I i : — ( II 

2 Ei Ej=il32 — j — ( e «(i) - e m>± + 4 — j ij ■ ( e tb) - e *w)^/' 



where 



D« =Dy + 3j(D x + D,)(e t y)-e tW ) 

= D y + I (E/=i D t(i)i' _ Ei»=i D t(i)i") ■ 
The first part of Eq. (|5"5|) is analogous to the expression (|49p . but with modificated DM vectors. It is important here 
that the combination (e t (j) —e t ^) varies from bond to bond with the use of the same symmetry transformations of the 
point group m3 as vector and, consequently, Dy. The vectors Dy are characterized by condition D x + D z = 0, so 
they have only two independent components. For example, the modificated DM vector of the bond (—2a;, \ ,\ — 2x) 
is 

D = (D x ,D y , D.) + ^(D X + D z )(e 2 «*) = fB^El,D v , ■ (60) 
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The second part of Eq. (|58j) contains the terms like 



aj_ • = a • b — (a • s)(b • s), 



(61) 



depending on the direction of s. The vectors a and b vary from bond to bond with the use of the symmetry 
transformations of the point group 23, and the averaging over 12 bonds in the unit cell, as well as for the case of the 
averaging over the sphere, gives 

12 , 9 

(62) 



1 12 12 
a ± ■ bj_) = a • b — ^((i?„a) ■ s)((R n b) • s) = a • b - -a ■ b = -a • b, 



12 

with R n being the rotation matrices of the group 23. Thereby, the mean contribution of the second part of Eq. ([5 
to the energy of a bond is 

D?) 2 , ,„ y/3D a + D a/rt „ „ , ,_ _ v, {D x + D z f 



2 3 (D x 

> I • 



(ea - ei ) 2 + x ' - (D x ,D y , D z ) ■ (e 2 - e x )} = -- 



3 l 32 J v 4 J 

and we can now rewrite the expression ()58[) for the energy as 

£ = \ E E{^(§ - 5) 2 + 6 « • (5 - §)} + (- 12 - 2 ' D ' D 

i j=l \ 



QJ 



J 



(63) 



(64) 



Let us introduce now smooth continuous vector fields Ci( r )i C2( r )j C3( r ) an d C4 ( r ) j connecting with Q by the equation 

(i = tm(*i)- ( 65 ) 
We can suppose that functions Q change slowly enough with coordinates and postulate the homogeneity condition 

d£ 2 ) n d(£ 3 ) n d(t)n _ d(n 



-, m, n = 1, 2, 3, 



(66) 



<9r m <9r m dr m dr m dr ri 

reflecting the fact that the divergence between sublattices is preserved over big distances. Then, we can write Eq. 
(|57p for 1st atom in a cell in the explicit form as 



Cx(x,x,x) = \ \ + x, \ - x) + C 2 (-x, -\ + x, \ - x) 

+C,z{\ ~ -x, \ + x) + £3(5 - x, -x, -\ + x) 
+Ca{\ + x, \ — x, —x) + Ci(— 5 +x,\—x, —x)^j 
= i ^C 2 (x,a;,x) + ( 3 (x,x,x) + ( 4 (x,x,x)j + | (^=1 r ij ' ^) CC^, a, a). 



(67) 



The last equation is correct not only for r = (x, x, x) but also for arbitrary r due to homogeneity condition. So, we 
can write Eq. (|67[) and analogous equations for other atoms in a cell, supposing that all Q are taken at the same 
point: 



(68) 



Ci =^(C 2 + C3 + C4) + i(E-=iryV)c, 
C2 =|(& + S + S) + |(E5=ir«-v)? > 
?3 = § (S + (2 + &) + | (E - =1 r 3j • V) <f, 

t S =i(Ci + C2 + c;)+i(E- =1 r4,-v)c*. 



The system (pjSf is degenerate on variables d, £2, C31 C4 an d satisfied, when 

6 6 



v 3=l 3'=1 

Let us allocate two parts from expression (1641) for the energy associated with a unit cell, 

J 4 6 J 4 6 

£ce«,e*cfc = J J2 Efe ( r i) - ^(r,)) 2 = J E E^« ~ 5 + ( F « ' W 



and 



i=l j=l i=l j'=l 



^46 ^46 

# ce «,zw = g E E 6 « ' (Ct(j)fa) - C<(rO) = 2 E E 6 « • (CtW) - Ci + (r« • V)C). 

i=l j=l i= 1 3=1 



With the use of Eq. (|69p we can rewrite these parts as 



2 

4 / / 6 \.\ T 46 



and 



BceH^cfc = E Ery • V C + ^ E E (( r « ' 

1=1 \ \J =1 / / i=1 3=1 



SceH^M = \ E E( r ^ • ^)(^ • cV g E ( E r « • v I IE ■ c 

»=1 3=1 *=1 \3=1 / = 1 



The second part of (|73[) vanishes, because 



Using 



6 

3=1 



=V3(l-8x) ei , 

3 = 1 



the first part of (|7^|) can be calculated, 



J &\ ?\ J ft o \2 d (n 9( n 



The second part of (|71?j) can be expressed as 

~ f; ^i?r„„((r 12 ) ; V fc a • ((ri 2 )„V m a = J(l - 4x + 16z 2 )|^|^, 

a— 1 

where summation is taken over all elements i? a of the point symmetry group to3, and we have used that 

24 

Epa pet 

and 

(r 12 ) 2 = i(f-4x+16x 2 ). 



OL=l 
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So, finally, 



Similarly, Eq. (|73|) can be rewritten as 



3J dC, n d( n 

E C ell,exch - JTTn-. (»Uj 

<± VI m Ul m 



E C ell,DM =^^mn( r 12)/(Dl 2 )nV fe Cm =4(r 12 -D 12 )-^ (81) 

Q— 1 

or 

E C eii,DM = —{D x ^ 2D y — D z )— — . (82) 

It should be noted that, whereas Eqs. (|72|) and (1731 contain parameter a;, the resulting expressions (l80l) and (|82|) 
do not depend on x. This looks like a coincidence, but it is not so. Indeed, the initial statement of the problem do 
not include atomic coordinates (e.g., see Eq. ([T])), and therefore the coordinates will not appear in results (at least 
until strong changes of the coordinates changing the topology of the net of nearest neighbors) . 

The substitution of Eqs. (J5DJ and flg2} into Eq. (J64J) gives 



with J and T> given by the following equations 



V = -4(r 12 • Dia) =D X - 2D y - D z . (85) 

The summation over all unit cells of the crystal has been replaced with the integration over the crystal volume. 
Let us suppose that the vector field £ is connected with a spin field s: 



s ss s + [C x s ], (86) 



where So is a constant vector, ( _L sq. Then, 



and 



So, we can rewrite Eq. ((55)) for the energy as 



ds n ds n ^ d( n 8( n 

dr m dr m dr m dr m 



(87) 



E = ! { jd & +m -v* s| ) * + N - J (- 12 - 2 (^) 2 ) ■ < 8<j) 

the first part of which coincides with the usually used energy in the phenomenological theory with the elastic constants 
J and V related to the the constants of microscopic interactions by Eqs. (|84[) and (l85|) . 
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V 2 (D x - 2D V - D z . 
k = —— = — I — , (90) 



The wave number for the helix structure can be now found using the phenomenological approach 

V _ 2 (D x - 2D 

which coincides with Eq. (1231) . 

Thus, we have passed from the discrete spin distribution of the microscopic theory to the continuous spin field s of 
the coarse grain theory. In order to complete our reasoning we can also consider the way back from the continuous 
approximation to dicrete spin distibution. 

Let s(r) be a continuous spin field from the phenomenological theory. Within a small domain of the crystal we can 
associate with s a vector field £ using the transformation inverse to 



C=[s xs]. (91) 
Here So is some mean value of s, which is constant over the domain. Then, 

^ = [soX — ]«[sx— ]. (92) 

Let us define continuous fields £i , C2 , C3 an d C4 f° r f° ur sublattices of magnetic atoms by formula 

c* = c+gE( r «- v K ) (93) 

obviously obeying Eq. Thus C is the average of these four fields, 

\& + & + & + &)- (94) 
Then, using Eqs. (1651) and (|5fc>[) we can express vector Q for individual atom trough £, 

Q = ?(r<) + ^ ^ + % (t)± (r a ) + ^(1 - 8x)(e t(l) ■ V)C>,), (95) 

where 

e t (i)x(ri) = e t (i) - s(r i )(s(r l ) • e tW ). (96) 
The substitution of Eq. (f95|) into the equation 

Si = s+[(,x s] (97) 

gives us the following expression for the spin of individual magnetic atom 

Si = S(r<) + ^. D * + D * [e t{i) x s(r,)] + ^(1 - 8x)(e t{i) ■ f )8(r t -). (98) 

This equation is valid for any smooth continuous field s(r). For one-dimensional helices, the first term in Eq. (|98l) . 
which is common for all sublattices, describes spin rotation in the plane perpendicular to the wave vector k, the 
second one defines a spin component parallel to k, and second and third ones taken together do a shift of the phase 
of the rotation for each sublattice. For instance, if one chooses s in the form of 

s(r) = ri! cos k • r + n 2 sin k • r, (99) 

then s, = s*W(r.j) has the form (pTTj) . 
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VI. DISCUSSION AND SUMMARY 

As it is shown above, in the approximation linear on (D/ J), only two components of DM vectors are significant 
for the considered magnetic structure of MnSi: D x — 2D y — D z and D x + D z . The third component, D x + D y — D z , 
appears in upper approximations. The pseudo-scalar constant T> used in the phenomenological theory is found to be 
proportional to the first component, D x — 2D y — D z . In fact, we could expect that only this combination of D-vector 
coordinates played a role in the simplest microscopic approximation. Surprisingly, it has been found that another 
component, D x +D Z , is also very important both because of its contribution into magnetic energy and its influence on 
spins orientations. This effect is very similar to "frozen" optical phonon displacements of atoms in crystals. Whereas 
the former component is responsible for the long scale chiral magnetic structure, the latter leads to a canting of spins 
within a unit cell and this canting could be in principle measured experimentally using diffraction methods t 23] . For 
instance, for the magnetic fields stronger than H Cl the residual splay of spins results in non-zero value of the neutron 
scattering amplitude for reflections of 00^, I — 2n + l type, forbidden both for nuclear scattering and for ferromagnetic 
scattering owing to 2i screw axes. 

An interesting result about the symmetry of the problem can be revealed considering the inverse structure, where 
magnetic atoms occupy positions ( ) rather than (x,x,x). The crystal thus obtained has the same space 

group, P2i3, but both the crystals are related as right- and left-hand structures, via the inversion transformation. 
It can be easily seen from the expression ((TJ that energy does not change over this transformation, if the spins in 
corresponding positions are equal and so do Dzyaloshinskii-Moriya vectors associated with corresponding bonds. This 
means that DM vectors are really pseudovectors. Also it is well understood H| that the inversion should change 
the chirality of the structure and, consequently, the sign of phenomenological constant T> and wave number k. However 
it is not so evident from Eq. that k changes its sign over inversion, because D is a pseudovector rather than 
a vector. We should rewrite the wave number in an invariant form, including an explicit dependence on the bond 
direction, in order to make transparent the change of chirality upon inversion: 

* = - 8(ri V 7 6l2) - (100) 

The pseudovector D12 does not change its sign under inversion but ri2 does. Notice that in this form the sign of k 
does not depend on our choice which of atoms is 1 or 2. 

It is interesting to compare our results with those ones obtained by Hopkinson and Kee [3] using mean-field 
and real-space rigid rotor minimization methods. According to this work, for both methods, a sizable spin helicity 
can be obtained only when the DM vectors lie parallel to the corresponding Mn-Mn bonds. Correspondingly, three 
components of DM vector are considered: the first one, parallel to the bond, the second one, perpendicular to 
the triangle of bonds, to which this bond is belonging, and the third one, lying in the plane of the triangle, but 
perpendicular to the bond. For instance, for the bond 1*12 = (—2 ) these three directions can be written 

in explicit form as (—2a;, h, \ — 2x), (1, 1,-1) and (1 — 2x, Ax — |, i + 2x), correspondingly, or, using x — 0.138, 
as (-0.276,0.5,0.224), (1,1,-1) and (0.724,0.052,0.776). The first component is responsible for the helix structure, 
and corresponding direction, (—0.276,0.5,0.224), is very close to direction (—1,2,1) for D x — 2D y — D z component 
found in our work (the angle between them is about 3.4 degrees). The second component, which does not cant spins 
[ijjl, coincides with our "neutral" component D x + D y — D z directed along the corresponding 3-fold axis. The third 
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component leads to spin canting and rather small spin spiraling, and its direction, (0.724,0.052,0.776), is very close 
to direction (1,0,1) for D x + D z component. Obviously, those numerical results almost coincide with our analytical 
solutions. An advantage of the analytical solution is that it does not rely on the exact bond directions and atomic 
coordinates; that is natural from the theoretical point of view because the atomic parameter x does not figure in Eq. 

m 

Moreover, different model considerations predict that the superexchange DM vectors are either perpendicular [24| or 
almost perpendicular 25H27| to the corresponding bonds. If this is also the case in MnSi, we could expect pronounced 



frustration of the helices. In any case, it is obvious that modern ab initio calculations of the Dyzaloshinskii-Moriya 
vectors would be very important for MnSi-type crystals. 

In conclusion, it is found which component of the Dyzaloshinskii-Moriya vector determines the long-range spiral 
ordering. The most interesting point is that the local structure of helices can be rather complicated owing to symmetry- 
driven frustrations. We hope that qualitatively all the considered features will appear also in more sophisticated models 
taking into account the itinerant nature of magnetism and thermal and quantum properties of spiral ordering in MnSi. 
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TABLE I: The bonds between neighboring Mn atoms in MnSi crystal and corresponding Dzyaloshinskii-Moriya vectors. The 
listed 12 bonds are connected by symmetry transformations of the point group 23. Other 12 bonds appear owing to permutation 
of i-th and j-th atoms {vji = — r^-, Dji = — ~Dij). 
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Figures 




FIG. 1: (a) Organization of a free helix with an arbitrary wave vector k (no magnetic field). The spins of each of four sublattices 
rotate around their own rotation directions n ro t (see Eq. (|30[0 . which differ slightly from that of k. In the considered model, 
the energy of the free helix does not depend on direction of k whereas experimentally the free helix is oriented either along 
[111] or [001] direction owing to the crystal field anisotropy. (b) Conical helix spin structure in an external magnetic field. The 
spins of each of four sublattices have a component directed along corresponding vector n rot ; its value is proportional to |H|/J? C 
for H| < H c and is equal to 1 for |H| > H c ; k- vector is directed along the magnetic field. 



22 



(a) 

mm • 

.Cv.(ioo) 



a 

fl 

o 

& 



(010) 

.V 4 V. (° 01 ) V. 



k || (001) 



360 
300 
4 240 



»U .V».....,J#«*S SMI a* g 



m • r •••• 



120 



- 60 



10 20 30 

r (unit cell parameters) 



40 



(c) 



c 

<u 

c 
o 
n. 

£ 
o 
o 
c 
'& 



. j^.(ioo) Tv.co, 

•A 3 (011) 'V, 



k|| (011) 

1.1) 




10 20 30 

r (unit cell parameters) 




10 20 30 

r (unit cell parameters) 



40 



(d) 




1 2 3 

r (unit cell parameters) 



FIG. 2: (a)-(c) Details of spin helix structures calculated with the use of Eq. (|3ip for the wave vectors k directed along (001), 
(111) and (011) crystal axes, correspondingly. The wave number |k| is chosen to be equal to 27r/40, which is close to the 
experimental value for MnSi crystal where the helix period is approximately 40 unit cells. The canting D-vector component is 
taken to be rather large, D x + D z — 0.4J, to emphasize nonmonotonic frustrated behavior of helices. Three spin components 
along directions m, 112 and ri3 = nk are presented by points with four colors corresponding to four sublattices. The black saw 
represents the angle <f) between spins in the successive ferromagnetic planes, tan0 = (s • ri2)/(s • ni). (d) An enlarged part of 
the saw restricted by a small box in (c); the line is drown to guide the eyes. 



